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Producing this graph required many calculations. The adjacent article 
has a method for storing the graph data for quick redraw. 



Solve Nonlinear Systems Of 
Equations In Seconds 

By Gary P. Laroff 

(This is the third in a series of articles on appUcations 
of the new Tektronix software package 4050A10 
Statistics Volume 4.) 

Equations that we are used to solving with pencil and 
paper are most often "linear" equations. Unfortunately 
many of the equations encountered in day to day 
experiences describe "nonlinear" phenomena, and are not 
so easily handled. The roots of a nonlinear equation f(x) = 
cannot in general be expressed in closed form. In order 
to solve nonlinear equations, we are thus obliged to use 
approximate methods. These methods, normally 
presented for a single nonlinear equation, can be 
generalized to systems of nonlinear equations. 

Consider the text book problem^: 

...The system of equations 

X = 1 + Wie'"" ' + 3x-) 
y = 0.5 + h" tan (e' -f y") 
(Continued on page 6.) 



Storing A Graph On Tape 
From The 4051 

By Dan Taylor 

Some 4051 applications require a large number of 
calculations to produce a graph (for example, a 3-D 
surface). It would be nice to be able to store the fmished 
graph on tape and get it quickly redrawn (on the screen or 
plotter) at a later time. This can be done in several 
different ways with the 4051 and its internal tape drive. 

The most straightforward method of storing a graph is 
to add new lines of code after each MOVE and DRAW to 
write (in binary) the X-Y coordinates on tape. But storing 
only the X-Y values is not sufficient because you wouldn't 
be able to tell the MOVEs from the DRAWs on the tape. 
A solution is to store a third number with each X-Y pair. 
The method I use is: 

1,X,Y for a MOVE to X,Y 

2,X,Y for a DRAW to X,Y 

This technique can then be extended to the other graphic 
commands: 

3,N,X1 XN,Y1 YN Arrav DRAW of N points 

4,X,Y RMOVE 

5,X,Y RDRAW 

6,N,XI XN,Y1 YN Array RDRAW of N points 

7 AXIS - no arguments 

8,A,B AXIS - 2 arguments 

9,A,B,C,D AXIS - 4 arguments 

10,\,AS PRINT the string AS of length N 

11, A PRINT the number A 

12,A,N,A$ PRINT USING AS: A (length AS^N) 

13,A ROTATE A 

14,A,B,C,D WINDOW A,B,C,D 

15,A,B,C,D VIEWPORT A,B,C,D 

16,A,B SCALE A,B 

A simple BASIC program with graphics is: 

100 VIEWPORT 30,100,20,80 
110 WINDOW -2,10,0,300 
120 AXIS 2,50 
130 MOVE 2,100 
140 DRAW 6,200 
(Continued on page 4.) 
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^ Editor's Note 



New 4051 Applications Library 
Catalog Included With This Issue 



The 405 1 Applications Library is now 8 months old. In 
that time we've collected 59 programs for business, 
engineering, scientific, educational and general utility 
applications. A number of these programs have come 
from members of our 4051 Applications Library 
members, and your support is greatly appreciated. The 
wide variety of 4051 applications demonstrated by these 
programs include making super 8 animated films, 
petroleum risk analysis, schematic drafting, economic 
analysis, educational demonstrations, as well as the 
complex mathematical analysis and graphing capability, 
which is the 405 Ts specialty. 



Do You Have A Friend Who 
Isn't A 4051 Applicationsr 
Library Member Yet? 

WeVe including an application card for membership in 
the 4501 Applications Library in this month's Tekniques. 
If you haven't sent us an application card yet, please do 
so, or give the card to a fellow 4051 user who would also 
like to receive Tekniques. 

In addition to receiving Tekniques, the 4051 
Applications Library newsletter, as a member of the 405 1 
Applications Library you are entitled to participate in the 
4051 software exchange program. Submit one 4051 
program to the 4051 Applications Library, and (if it is 
accepted) receive three programs of your choice from the 
librarv for free. 



With this issue of Tekniques is a 12 page catalog of the 
programs collected to date by the 4051 Applications 
Library. This catalog includes reprints of the program 
abstracts that have appeared in the past issues of 
Tekniques. an index of the programs, and ordering 
information. All programs are available in our software 
exchange program. The catalog will be updated 
periodically. Additional copies of the catalog, are 
available from your local Tektronix Sales Engineer, or by 
writing Tektronix, Inc. P.O. Box 500, Group 451, 
Beaverton, Oregon 97077 c/o Applications Librarv, (60- 
369) 

Besides developing a library of software for the 405 1 , 
we also see the 4051 Applications Library as our most 
effective means of interfacing with you the 4051 user. 
Submitting a program to the 4051 Applications Library 
not only allows you to share your efforts and knowledge 
of the 4051 with other 4051 users, but it also gives us a 
better idea of how we can meet your needs in the future. 
Even if you don't have a program to submit, we are always 
interested in hearing how you are putting the 4051 to 
work for you. 
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An Unusual Use of MOVE and DRAW 

By Dan Taylor 

In normal usage, the key word MOVE transforms 
coordinates from user units (window units) to graphic 
display units (GDUs). If the GDU location is off screen, 
nothing happens. If the GDU location is on screen, then a 
move (with the beam ofQ is done to the point, even if it is 
outside the viewport. 

DRAW also transforms user units to GDUs. A vector 
is then drawn toward the point. If the point is outside the 
viewport, the vector is stopped (clipped) at the viewport 
boundarv. 



to: 



The simplest code for a draw and a move is equivalent 

DRAW (5)32,20: 
MOVE@32,21: 

Where 32 is the address of the screen, 20 tells the screen to 
move with the beam on; and 21 tells the screen to move 
with the beam off. 

The screen (or plotter) doesn't know what keyword 
(MOVE or DRAW) was called, only whether to move 
with the beam on or off. You can thus call a MOVEl^ 
@32,20: to draw a vector that goes outside the viewport.^ 
Similarly, you can call DRAW @32,21: if you want a 
move (with the beam ofO clipped at the viewport 
boundarv. 



^ 



GOTO And GOSUB Extend 
Branching Capability 

By Ken Cramer 

If you are an ex-FORTR AN programmer, some of the 
BASIC constructs in the 4051 may seem limited, par- 
ticularly the branching associated with the true and false 
consequences of IF statements. The computed GOTO 
(GOTO n OF 500,600,700...), however, in 4051 BASIC 
allows extended branching capabilities similar to FOR- 
TRAN. 

The FORTRAN IF syntax allows branching to three 
different statement numbers as a consequence of a minus, 
zero, or positive result. When using the 4051, the 
FORTRAN statement: 

IF(A) 100,200,300 

can be rewritten in 4051 BASIC as: 

299 GO TO £GH<A)+2 OF 588,608,798 

599 REM BRANCH GOES TO HERE IF ft IS HEGATIUE 

519 REM . 

539 REM CONTENTS OF ROUTINE 

549 REM . 

558 GO TO 1999 

699 REM BRANCH GOES HERE IF A IS EXACTLY ZERO 

618 REM . 

638 REM CONTENTS OF ROUTINE 

659 GO TO 1989 

789 REM BRANCH GOES HERE IF A IS POSITIUE 

718 REM . 

738 REM CONTENTS OF ROUTINE 

749 REM . 

798 GO TO 1898 

If A is an expression or A is a computed value, the chances 
of A being exactly zero are pretty small. To bring the 
FUZZ factor into the evaluation of A for a negative, zero, 
or positive result, use the following statement: 



When you want to execute a subroutine on a true or 
false condition, a dummy RETURN statement can be 
used to effect the same result. For example, the following 
FORTRAN operation: 

1F(A<0) CALL SUB 1 

can be accomplished in 4051 BASIC as 



A''9>-^1 OF 599,698 
THE SUBROUTINE RETURNS HERE 



299 GOSUB 

218 REM 

229 REM 

239 REM . 

599 REM THE SUBROUTINE STARTS HERE 

519 REM . 

529 REM . 

599 REM 

595 REM 



699 RETURN 



THE FOLLOWING RETURN STATEMENT IS ALWAYS EX 
REGARDLESS OF A's UALUE 



Easy Cyclical Counters 



By LeRoy Nollette 



It is often necessary to have a counter in your program 
that counts up to a certain number, then starts over. A day 
of the week counter, which counts up to 7 then starts over, 
is one example. The simplest way to write such a counter 
is: 

190 W=0 

110 W=Nt(W<7)^l 



299 GO TO SGHCA)*<A«9)+2 OF 588,698,789 

588 REM BRANCH GOES HERE IF A IS HEGATIUE 

518 REH . 

539 REM CONTENTS OF ROUTINE 

549 REM . 

599 GO TO 1899 

699 REM BRANCH GOES HERE IF A IS WITHIN FUZZ OF ZERO 

619 REM . 

639 REM CONTENTS OF ROUTINE 

699 GO TO 1898 

789 REN BRANCH GOES HERE IF A IS POSITIUE 

719 REM . 

738 REM CONTENTS OF ROUTINE 

748 REM . 

799 GO TO 1899 



160 GO TO 110 

When W in line 110 is less than 7, it is increased by 1. 
When W is equal to 7 it becomes 1. This counter can be 
written to count to other numbers merely by changing the 
7 to the desired upper limit. 



J 



4051 BASIC expands the computed GOTO capability 
with the computed GOSUB statement. GOSUB is not 
available in FORTRAN. The FORTRAN statements: 

IF(A=1) CALLSUBl 
IF(A=2)CALL SUB2 
IF(A=3)CALL SUB3 

can be reduced in 4051 BASIC to one line of code: 



299 GOSUB SGN<A>»<A«9)+2 OF 588,688,799 

219 REM EXECUTION OF ANY OF THE SUBROUTINES RETURNS HERE 

229 REM 

599 REM BRANCH GOES TO HERE IF A IS HEGATIUE 

518 REM . 

538 REM CONTENTS OF ROUTINE 

548 REM . 

599 RETURN 

688 REM BRANCH GOES HERE IF A IS WITHIN FUZZ OF ZERO 

619 REM . 

639 REM CONTENTS OF ROUTINE 

698 RETURN 

789 REM BRANCH GOES HERE IF A IS POSITIUE 

719 REM . 

738 REM CONTENTS OF ROUTINE 

748 REM . 

799 RETURN 



GPIB Manual Now Available 



The GPIB Application Support manual is now 
available through your Tektronix Sales Engineer. The 
part number to order is 070-2307-00 and the price is 
$20.00. This 250 page manual discusses the IEEE 
standard, the signals on the bus, data transmission and 
the GPIB handshake. Major sections of the manual are 
dedicated to program examples using Tektronix, Dana, 
Fluke, Wavetek and Hewlett-Packard equipment. 

This manual will be useful to everyone interested in 
using the 4051 GPIB to interface to IEEE-488 equipment 
built by other companies. 



storing A Graph ■ ■ - ■ from page 1 

To store this graph on tape you would add the following 
lines to the program: 

95 FIND J' 

100 VIEWPORT 30,100,20,80 

105 WRITE 15,30,100,20,80 

110 WINDOW -2, 10,0,300 

115 WRITE 14,-2,10,0,300 

120 AXIS 2,50 

125 WRITE 8,20,50 

130 MOVE 2,100 

135 WRITE 1,2,100 

140 DRAW 6,200 

145 WRITE 2,6,100 

147 CLOSE 

Note that only the code numbers for the graphic 
commands and coordinates are stored on the tape. The 
calculations required to create these coordinates do not 
have to be performed again. The graph, or 3D surface can 
thus be redrawn very quickly from this information. 

The following program can then be used to read such a 
tape file and reproduce the graph: 



lit IS out 



32> 



leee PPIHT '•L''«-P'*a»« » ^raoh unich is store'J 
leia PPIHT "JEnter the Uoc file « 'I 

t fl'^i^ I MP tJT J 

1939 IF J<VIHT<J^ OR J< 1 OR J>255 THEN 1919 
19-49 PRINT -JEnter th« device nuwtoer for jra 

[III iF''D0<>?HT<D8^ OR 09< 1 OR D9>33 THEN 1848 

1978 FIHD J 

1988 OH EOF <8) THEN 4889 

1189 REM GET THE NEXT GRAPHIC ELEMENT 

II28 Go'^TO*"oF 2189,2290,2398.2499,2599.2699.2799,2899.2999.3999 

IU9 GO 70 K-ie OF 3198.3299,3393,3499.3399.3699 

1149 PRINT 'LError: Th;a fraoh * '- 

1159 END 

2998 REH HOUE: 

2198 READ f33:X9,Y8 

2118 NOUE eDe:xe,Y9 

2128 GO TO 1199 

2199 REH dram: 
2299 READ •33:X9,Ya 
2218 DRAW fDe:X8,Y8 
2228 GO TO 1188 
2298 REH MATRIX DRAM: 
2398 READ f33:H8 
2318 DIM Xl<M9),Yl<Ne> 
2328 READ fSS'.XliYl 
2338 DRAW tD0:Xl,Yi 
2348 DELETE XI, Yl 
2358 GO TO 1189 
2398 REM RHOUE: ^ 
2498 READ tSStxa.Yd 

2418 RMOWE fDe:xe,Y9 

2429 GO TO 1100 

24'?9 REM ROPAW: 

2599 REhD l?33:XO.Y0 

2519 RDRt^W 9D9:X9, V9 

2529 GO TO 1199 

2599 REM MATRIX RDRAW: 

2699 READ 833: N8 

2619 DIM X1<H8),Y1<H8> 

2629 READ t33:Xl,Yl 

2639 RDRAW »D8:X1,YI 

2649 DELETE XI, Yl 

2659 GO TO 1198 ^ ^^ 

2698 REH AXIS — HO ARGUMENTS: 

2?99 AXIS 9De: 

2719 GO TO 1188 

2799 REM AXIS — 2 ARGUMENTS: 

2899 READ »33:xe,Y8 

2818 AXIS »De:X8,Y8 

2828 GO TO 1188 

2998 REH AXIS — 4 ARGUMENTS: 

2988 READ f 33:X8, YB, X2, Y2 

2918 AXIS fD8:Xe,Y9,X2,Y2 

2928 GO TO 1188 

2998 REH PRINT A«l : 

3988 READ ff33:H8 

3819 DIM A»<M8*1> 

3828 READ »33:AI 

3838 PRINT •De:AS| 

3848 DELETE Af 

3838 GO TO 1108 

3898 REM PRINT Al : 

3188 READ t33:X8 

3118 PRINT fDe:X8l 

3128 GO TO 1188 

3198 REM PRINT USING Ai:A : 



The above technique is very general. There are a 
number of more restrictive methods of storing the graph 
on tape that can greatly increase the speed of the 
redrawing. One factor that slows redrawing is that the 
commands MOVE and DRAW X,Y must transform the 



X,Y values from user defmed units into graphic display 
units (GDUs) and perform any clipping required at the 
VIEWPORT boundaries. . 

If your graph is produced with just the WINDOW, 
VIEWPORT, MOVE, DRAW and PRINT commands, 
and no clipping occurs, it is fairly simple to speed up the 
redrawing, by converting user coordinates to GDUs 
before they are stored on the tape. 

If you use WINDOW W1,W2,W3,W4 and VIEW- 
PORT PI, P2,P3,P4, then you can construct the follow- 
ing relationships to convert user defined coordinates to 
GDUs: 

a = P2-Pl 



W2-W1 

b = PI - a*Wl 

c ::: P4 - P3 
W4 - W3 

d = P3 -c*W3 

X and Y, in GDUs, is thus equal to: 

X = a * X + b 

Y = c * y + d where x and y are in user defmed 
units. 

Note that MOVE x,yis now the same^ as PRINT @32,21jr 
a * X + b, c * y + d, and DRAW x,y is the same as PRINT^ 
@32,20: a ♦ X + b, c * y + d. 

As with the first method, the graphic data can be 
written on tape in triplets, in the form: 



20,X,Y 

21,X,Y 

1,N,A$ 

2,A 

3,A,N,A$ 



draw to X,Y (X,YareinGDU's, not 

movetoX,Y user units) 

print AS of length N 

print the number A 

print USING AS: AS is of length N 



In writing a program to create this tape file, you use the 
relationships for the conversion factors a, b, c and d to 
convert coordinates from user defined units to GDUs. 
For example: 



3209 READ (?33:N9,N0 

3219 DIM A«vH0*r^ 

3220 READ 833: Al 

3239 PRINT 9D0; USING Ai:x9 

3249 DELETE A« 

3250 GO TO 1109 
3299 REM ROTATE: 

3398 READ •33:X9 
3310 ROTATE X9 
3329 GO TO 1100 

3399 REM UIMPOW: 

3499 READ f33:X8,Y8,X2,Y2 
3419 UlNDOU X9,Y0,X2,Y2 
3428 GO TO 1189 
3490 REM UIEHPORT: 
3398 READ #33: X8, Y8, X2. Y2 
3518 VIEWPORT X8,Y0.X2.Y2 
3528 GO TO 1188 
3598 REM SCALE: 
3688 READ t33:X8,Y8 
3618 SCALE Xe,Ya 

3998 REH^McTURC FINISHED - MOUE CURSOR TO UPPER RIGHT: 
4888 INPUT ffD8:X8,Y8 
4818 PRINT tD8,2t:X8iY8 
4828 HOME 
48 18 END 



( 



The program to read the resulting file from tape does all 



") 



) 



J 



the graphics with PRINT commands, using the 
parameters 20 and 2 1 to determine if the print is a move or 
a draw (see the 405 I Graphic System Reference Manual, 
pages 9-27 and 9-4 1 , for the use of secondary addresses 20 
and 21). 



U0e PEflO 933:5 

nia GO ro s of 3090.3108,3209 

1129 PEM HOWE OR DRAW VGDU LEUEL^ 

1139 READ »33:>:0.V9 

1149 PRINT er'0.s:xe. Y0 

1159 GO TO 1 199 



It is several times faster than the first method. If your file 
doesn't have any PRINT commands this simplifies to: 



1199 PEhC) i?33:s.xe,vi 

I 1 19 PRINT ?09,S:X9, 
1 129 GO TO 1199 



The last technique is another two to ten times faster 
than the second method. It makes use of the fact that a 
graph normally consists of a move followed by a number 
of draws. The structure on tape is: 

K,Z (forK^ 2 and even). K is the length of the linear 
array Z. And Z consists of alternating X and Y 
values - move to the first X,Y pair in Z and then 
draw to the rest. (X,Y are in GDU's) 

-1,N,AS Print, as previously 

-2, A Print, as previously 

-3,A,N,A$ Print, as previously 

A program to read such a tape file and produce a graph is: 



1199 READ ^33:). 




I I 19 GO TO -H OF 3990 


310 


1129 REM MATRIX "CUP^ 


'E- 


1139 DELETE 2 




I 149 DIM Z<<) 




1159 READ g33:Z 




1169 PRINT ?D0.2l:2a 


, z< 


ll?9 IF K-2 THEN 1198 




1188 PRINT »D9.29:2 




1199 GO TO 1188 





When there are an average of 10 draws for every move 
this method is 7 times faster than the second method. The 
main difficulty with this last method is that your program 
to create the appropriate tape file becomes more com- 
plicated. (Because you have to know how many draws will 
follow before any can be put on tape.) 



'The size that the tape should be marked depends upon how complicated 
your graph is. Each number in a binary file requires 10 bytes. Each string 
requires 3+ the number of characters in the string. A normal graph might 
have 100 draws (3 x 100 x 10 = 3000 bytes), therefore, using 10,000 byte 
files would be adequate for most graphs. (Binary files are used because 
they are faster.) 



-Except for clipping i.e. this is true as long as Wl ^ x^ W2 and 
W3^^ y ^ W4. 



More Information on 
Computer Animation Using 
the 4051 

By Will Gallant 

Last month's article by Mr. Paul Doherty of Oakland 
University in Rochester, Michigan on computer anima- 
tion with the 405 1 (see Tekniques Vol. 1 No. 5) presented 
some useful techniques for making super-8 animated 
films with the 4051. Several of us have been experimen- 
ting with computer animation and our after hours 
recreation in cinema photography has produced some 
notes we'd like to share with you. Like Mr. Doherty we've 
discovered that with some method of automatically 
triggering the camera single frame and some thoughtful 
programming, you can make animated films quickly and 
efficiently with a minimum of operator interaction. 

A camera we use is the Minolta XL400. It has a single 
frame electronic shutter and a good /1. 2 lens. 

We use a different method of triggering the electronic 
shutter than suggested in the last issue of Tekniques. Pin 
1 1 of the hard copy jack (J41) on the 4051 provides a 1/25 
of a second trigger pulse whenever a COPY command is 
executed. The single frame input on the Minolta XL400 
triggers when its input drops below 300 ohms. Fig. I 
shows the interfacing required to trigger your camera 
from the hard copy jack. 




>-[ 




J41 Hdrdcopv Jack 

(Pins n'dnd 15) 



Minolta XL400. 
using ASA160 
or Faster Fr.'^. 



Fig. 1. Interfacing required to trigger camera's electronic 
single frame from hard copy jack. 

Also as Mr. Doherty indicated, the resolution and 
brightness of the 4051 pushes low cost super 8 cameras 
and film to the limit. If you can tolerate less detail, movies 
made with Kodak ASA 160 color film are adequate for 
projection in a very dark room. This film is readily 
available and adequate for most filming of graphic 
objects, however, text is very difficult to pick up clearly. 
GAF makes an AS A500 black and w hite film for super 8, 
but again you can expect some loss of detail in the 
bargain. 

For maximum brightness of the movie, macro- 
focusing lenses are required, but they are hard to focus 
because of the small depth of field. One trick we learned is 
to hold a ruler perpendicular to the glass on the 4051 
screen. The reflection of the ruler on the screen allows you 
to read the distance from the glass to the phosphor 
(Continued on page 8.) 



Sol ve NonrmearEqu, ,.■■■. from page 1 

can, for small values of h, be solved by iteration. Put x,. 
- 1 ,' yo = -0.5. These values are put into the right-hand 
side, where upon a new approximation (Xi, yi) is 
obtained. The process is repeated until the changes in x 
and y are sufficiently small... solve the system for h - 0, 
0.01, and 0.10. 

Program 4 of Statistics Vol. 4, Nonlinear Systems Of 
Equations, can solve such problems, and is adaptable to 
the user's environment because it has been segmented into 
three parts: 

1. A subroutine to solve the system of nonUnear 
equations, 

2. a main program to use the subroutine, and 

3. another subroutine, user written, that defines the 
system of equations. A sample subroutine is 
included in the initial program printout. 

Method of Calculation 

PowelFs algorithm, used in Program 4, converges 
rapidlv with nearly global convergence. The algorithm 
combines features of the Newton-Raphson, the 
Levenberg-Marquardt, and the steepest descent methods. 

Nonlinear Svstems of Equations (NLSE) requires N 
equations of the form f,(Xi, X: ..., Xn) = 0. N represents 
the number of variables, and therefore the number of 

equations.]- 1,2 N; and N ^ 1. In vector form. X==(x,, 

..., xs), and F(X) - (f,(X),f:(X),..., fslX)). If N=1,and 
Xk ==^ value of x at kth iteration, Newton s method 
would be: 

x.H =^ X. - f(xo = X. - {r(x.)r * f(XO. 
f'(xo 

if N>2. Xkh = X, - J"'(X0*F(X0. 
J( Xk ) is the Jacobian matrix of F, evaluated to Xk . 

Variable D0 is used for approximating partial 
derivatives of the functions 



9x, 



(f,(X-}-D0 * e.) -f,(X))/D0 



J„ = partial of f. with respect to x,, evaluated at 



X,= 



„ dl 



9x, 



Jij is approximated numerically by: 
(f,(X,+D0*e,) -f.(XO)/D0 

where e, is the normalized jth coordinate vector = a 
vector except a I in the jth component. 

Variable E represents required accuracy, and the 
algorithm has found a solution, X, when: 

N 

1 (f,(X))^^E. 

P1 



The algorithm takes the sum of the squares of the f (X) 
values. When this sum is less than a specified E value, the 
program stops. 

Using The Program 

The user-written subroutine, starting at line number 
1 100, must define the system of equations and must end 
with a RETURN statement. Using the text book example 

above, 

1. Rearrange equations to get all terms on the right- 
hand side: 

= 1+ h'(e^''^ + 3x') - X 
== 0.5 + h'tan(e' + y") ^ y 

2. Set the first equation equal to F(l), the second to 
F(2),....and the nth the F(n), 

3. Replace the variables with x(i): 

X(1)-x 

X(2) - y 

The equations now look like: 

FU)-UHt2<CEXP(X(2>tSQRX(l)>*3tX<l>t2)-X<l) 
F(2)=e.5*Ht2lTAH<.EXP(X(l>>>XC2)t2)-X(2> 

Code these equations into the user-written subroutine 
starting at line 1 100, as requested by the main program. 
Replac^e H in these equations with 0, 0.1, or 0.01 as you 
key in the equations. 

For H - 0, the user-written subroutine looks like this: 



c 



II29 Pe! i"""-isAMp[ruiErwRITTEHluBRrTO EUrtLUnTE f AT :<1- 
1140 F. 1 . = Udt2t^EyP^X(2-'tS0RCM«;i>y'*3tXO^T2)-X* 1^ 



♦erstTftHiE'^P'; 



lire REN I- tEHD OF SAMPLE USER HRITTEH SUBROUTIHEt- 

1189 REH X " ""' 

1133 REM 



For H :=0.l: 



II29 re" r-------«S;^NPLE"iER WRITTEN SUBR. TO EUALU^TE F PiT X.-^ 

1149 FO n = 1 ♦e . 1 T2t (-EXPt X ( 2 > tSQR' X ^ i > >] ♦StX ^ t2 ' -XC 1 > 
1138 F.-2-8.5*9.ir2«THHfEXPtXCl'-'fX>2./T2'-X<2^ . 



lira SH^t- -tEHP OF SAMPLE USER URITTEH SueROUTINEt- 

1189 REM I '■ -" 

1198 REM 



And for H = 0.01: 



1189 REM ^_ , 



p1 pin i tSAMPLE USER WRITTEN SUBP. TO EVALUATE F AT X« » 

U9 F(l^'l*9.91T2»CEXP<X(2-.»SQRiX<l)>)*3tX<lU2^-XO. 
1159 F.-25=9.5*9.9lT2»TAN(EXPU(P^*X(2)T2-'-X(2^ 

1178 rIm^I- tEHP OF SAMPLE USER HRITTEH SUeROUTIHE* "^j 

1189 REM I— --.r.- - — - 

1199 REM 

To solve for H = O.Ol, key in RUN after the user- 
written subroutine, and enter the initial parameters for 
X(l) and X(2), 1 and 0.5 respectively. Then enter the 
desired accuracy. In this case the iterations are to be 
broken off when the sum of the squares of the function 
values are less than 10-^. 

NONLINEAR SYSTEMS OF EOUATIOHS 

ENTER N. TN£ NUMBER OF EQUATIONS - 2 

ENTER THE u£CTOR X, THE INITIAL UALUE FOR ITERATIONS: 

XI • 2 
X2 - I 

ENTER E, THE ACCURACY REQUIRED IN THE SOLUTION vEC, 1E-19> ■ lE-B 

ENTER M, THE "OELTA" FOR THE DERIMATIUES vEG, lE-3> • lE-3 

ENTER THE MAX. • OF ITERATIONS TO BE ALLOWED • 198 

INTER M. AH UPPER BOUND ON THE DISTANCE TO A SOLUTION <EG. IM) - 189 

DO YOU WANT X t F PRINTED AT EACH ITERATION STEP: YES 



C 



The time to solve the equations is 5 seconds, including 
the time to print out the intermediate results: 



1 



^ 



AFTER 


I 


EUALUTIQHS OF F: 
J '.« J) 

1 1 

2 9.5 


f J> 

4.64972127IE-4 

-l,73867l549E-5 


AFTER 


2 


E'lALUTIOHS OF F: 

1 i.asi 

2 9.5 


FU> 

-5.3449636E-4 

-1.722654876E-5 


AFTER 


3 


EVALUTIOHS OF F: 
J A(J> 

1 1 

2 9.581 


F*" J) 

4.659379817E-4 

-0.98191749356362 


FIHfiL ! 


SOLUTION, AFTER 4 EUALUTIOHS 
J :<(J> 

1 1.99946516767 

2 e. 499982621816 


OF F is: 

F(J) 

-7.289457926E-11 

3.694822226E-13 



.> 



Even with a choice of initial values of X(l) = 2 and X(2) = 
1, the total time ofthe calculation is only 5 seconds. For H 
- 0, the computation time is well under I second. 

In the next few months, look for part 4 on applications of 
Fourier Regressions to nonlinear function analysis. 



'Numerical Methods by Germund Dahlquist and Ake Biorck. Prentice- 
Hall, Englewood Cliffs, 1974, page 253. 



4051 Applications Library 
Program Abstracts 



Documentation and program listings of these 
programs may be ordered for S15.00 each. Programs will 
be put on tape for an additional S2. 00 handling charge per 
program and a S26.00 charge for the tape cartridge. (The 
program material contained herein is supplied without 
warranty or representation of any kind. Tektronix, Inc. 
assumes no responsiblity and shall have no liability, 
consequential or otherwise, of any kind arising from the 
use of this program material or any part thereof.) 

Please use the Applications Library Order Form. 
Order forms are included in the Membership Packet and 
are available from your local Tektronix Sales Engineer. 

• SCIENTIFIC 

ABSTRACT NUMBER: 51/00-5401/0 

Title: Graphics Demonstrations For Astronomy and 
Phsyics 

Author: Dr. R.J. Reimann, Physics Department, 

Boise State University 
Memory Requirement: 8K 
Peripherals: None 

This is a program series, under the control of a 
directory, for the presentation of introductory astronomy 
and physics concepts to students with no computer 
experience. The series includes: 

Sidereal Time— an aid to understanding. 

Solar Time— a view of circumpolar constellations at 

any time. 
Planetary motion — Circular orbits for the inner 

planets are plotted using Kepler's third law. 
Binary Stars — circular orbits for binary systems are 
plotted using Kepler's third law. 



Satellite Orbit— earth satellite orbits are plotted using 

a 2-step iterative approximation of Newton's law 

of universal gravitation. 
Stellar Magnitudes — relative apparent magnitudes of 

stars are drawn assuming that brightness is 

proportional to area. 
Phaser — a demonstration of the reference circle as 

related to simple harmonic motion: 
Vector Combinations — addition of any number of 

vectors by the head-to-tail technique. 
Trajectory With Drag — a quick 2-step iteration 

assuming drag is proportional to v^. 
Linear Least Squares Fit With Plot — best fit to a 

straight line is determined, 

ABSTRACT NUMBER: 51/00-5701/0 

Title: Kaplan-Meier Survival Table Compu-Plotter 

Author: Paul W. Baim, Division of Biostatistics, Univer- 
sity of Miami School of Medicine 
Memory Requirement: 24K (Option 21) 
Peripherals: 4631 Hard Copy Unit (Optional) 

The survival table is a means of measuring the response 
or nonresponse of subjects over time. This response may 
mean improvement of patients receiving a new drug, or 
failure of light bulbs using a new type of filament: both 
lend themselves equally well to this type of study. 

The program accepts raw data in the form of dates and 
subject status. The dates correspond to the beginning and 
end of the interval during which a particular subject is 
under study. The subject status indicates whether the 
subject is still responding at the end of the interval, or has 
stopped responding; for instance, whether a battery is still 
producing current after six days time, or has gone dead. 
The program then converts the raw data into a table of 
three columns: the interval, the cumulative survival, and 
the standard error. The interval is the same as that 
mentioned earlier, the cumulative survival represents the 
percentage of subjects surviving (i.e. responding) for this 
particular interval or longer, and the standard error 
establishes the reliability of the cumulative survival 
figure. Next, the program will, on demand, graph the 
calculated values, with the time inter\'als on the X axis 
and the cumulative survival percentages on the Y axis. 
The standard errors are shown by vertical lines. The 
purpose of the graph is to show clearly, how subjects are 
responding. Finally, the program will, on demand, list the 
data as it was entered. 

• MISCELLANEOUS 

ABSTRACT NUMBER: 51/00-6002/0 

Title: Print Mail Addresses and Form Letters 

Author: Nick Fkiaras 

Memory Requirements: 24K (but 32K is recommended) 

Peripherals: 4641 

This package actually contains two programs. 
Program One 

This program allows the user to enter the text of a form 
letter and save it on file 3 of the same tape. The letter can 
have several pages, which are numbered automatically. 
The text may be edited by adding, inserting, deleting and 
replacing lines. Once the letter has been composed and 



edited, it can be printed using the second program in this 
package. 

The user may have the same salutation printed 
automatically for every addressee (entered and saved with 
the text of the letter) or he may enter it himself for each 
addressee as the letters are being printed. 

Program Two 

This program allows the user to enter any number of 
mail addresses, save them on tape automatically, then 
recall them later and have them printed on "Pin Feed 
Labels" with a Tektronix 464 1 Printer. Each mail address 
can have up to five lines. Each line can be up to 32 
characters long. 

The user can edit his "mail address" files by deleting 
any number of mail addresses from them and by 
changing, replacing, inserting or deleting lines from each 
mail address. A permanent record of ail files may be made 
by printing them on the 4641 printer using regular printer 
paper (132 columns). 

Finally, the program will print form letters using the 
mail addresses and the text entered and saved by program 
one. More than one "mail address" file may be used. 

• GRAPHICS 

ABSTRACT NUMBER: 51/00-9501/0 

Title: Dash, Dot, Dash-Dot Routine 

Author: Nick Fkiaras 
Memory Requirement: 8K 
Peripherals: 4662 (Optional) 

The program draws a dotted, dashed, or dot-dashed 
line between any two points XI, Yl and X2,Y2 regardless 
of the window and viewport used. An example illustrating 
how to use the routine is included. 

ABSTRACT NUMBER: 51/00-9513/0 

Title: Presentation Aids, Release 2 

Author: Will Gallant 

Memory Requirement: 8K 

Peripherals: 4662 Interactive Digital Plotter 

This program is designed to aid in making 8 1/2 x 1 1 
inch overhead transparencies. Text may be input to the 
4051 and drawn on the 4662 Interactive Digital Plotter 
with full control over position and letter size. Special 
functions such as left and right justification, lines, boxes, 
and diamonds are also included. 

A pie chart routine is also included in this program. 
Any number of segments may be assigned text and a 
value. The annotation is printed radially on the center line 
of its respective sector of the pie. the program will 
compute and annotate each segment in percentage form if 
desired. Further, any one segment of the chart can be 
exploded a specified distance from the pie for graphic 
emphasis. 

ABSTRACT NUMBER: 51/00-9514/00 

Title; Pie Chart 

Author: P.C. Holman, University of Wisconsin — Stevens 

Point 
Memory Requirement: 16K 
Peripherals: 4631 



Data to be presented as a pie chart is entered on. the 
4051 keyboard. The data is then presented as a table of 
original data, percent conversions, and degrees of a circle 
conversions. Data is automatically presented as a pie 
chart in one of two sizes, in one of seven screen locations. 
The program is tutorial from tape (ASCII). 



Computer Animation , , . , , fro m page 5 

(typically 3/4 inch). This can then be added to the 
distance from the glass to the lens. 

If you can tolerate dimmer movies, it is much easier to 
place the camera back about 5 feet from the screen and 
zoom in. Your depth of field in this case increases to about 
M/2 inches. For maximum brightness of the finished 
movie, the film should always beshotina very darkroom 
or construct a tube between the 4051 and the camera. 

With the availability of the automaUc shutter control, 
you should be able to write programs to create the graphic 
images to be filmed without the need to have an operator 
on hand. The following is a simple program that suggests 
the general approach. Note that after the COPY com- 
mand in line 240, the screen is erased eight times. This 
erase loop is necessary because of the automatic origin 
shift in the 4051. It is recommended that this feature not 
be disabled. The total time to complete this erase loop is 
only 6 seconds, which doesn't add much to your per frame 
shooting time. 

Here are a few methods of speeding up filming. When 
you make three dimensional transformations and 
rotations, the opUonal Matrix ROM pack may speed up 
calculation time as much as 3 to 10 times. You might also 
do your calculations before filming and store the graphics 
on tape. (See the article on page 1 of this issue of 
Tekniques titled "Storing Graphics From the 4051 On 
Tape",) 

Another technique found to speed up filming is to take 
two frames of each graphic display, and then use larger 
increments of change between each transformation or 
rotation of the graphic image. You get the same rate of 
change of the image, but it takes half as long to film. 
Remember to wait a half second or so between COPY 
commands to allow your shutter to cycle. The following 
table gives you the filming and projection times for 
various intervals between frames. 

Filming And Projection Times 



Filming interval Filming time 
between successive required to 
frames in seconds, produce 10 
seconds pro- 
jection time 
at 18 fps. 



Time required to 
expose 50 feet of 
film. 



.5 


1.5 min. 




1 


3 min. 


1 hr 


2 


6 min. 


2 hr 


4 


12 min. 


4hr 


8 


24 min. 


8 hr 


15 


45 min. 


15 hr 


30 


1 hr. 30 min. 


30 hr 


60 


3 hr. 


60 hr. 



c 



